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We present cosmological perturbation theory in neutrino probe interacting dark-energy 
models, and calculate cosmic microwave background anisotropies and matter power spec- 
trum. In these models, the evolution of the mass of neutrinos is determined by the 
quintessence scalar field, which is responsible for the cosmic acceleration today. We con- 
sider several types of scalar field potentials and put constraints on the coupling parameter 
between neutrinos and dark energy. Assuming the flatness of the universe, the constraint 
we can derive from the current observation is J^m^ < 0.87ey at the 95 % confidence 
level for the sum over three species of neutrinos. 
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1. Introduction 

After SNIeAj and WMApEI observations during last decade, the discovery of the 
accelerated expansion of the universe is a major challenge of particle physics and 
cosmology. There are currently three candidates for the Dark-Energy which derives 
this accelerated expansion: 

• a non-zero cosmological constantP, 

• a dynamical cosmological constant (Quintessence scalar fieldj^, 

• modifications of Einstein Theory of Gravitj^ 

In this paper, we review shortly the main idea of three possible candidates and 
their cosmological phenomena. Specially we consider the interacting mechanism be- 
tween dark-energy with a hot dark-matter (neutrinos). Within neutrinos probe in- 
teracting dark-energy scenaricP^ we calculate Cosmic Microwave Background(CMB) 
radiation and Large Scale Structure(LSS) within cosmological perturbation theory. 
The evolution of the mass of neutrinos is determined by the quintessence scalar 
filed, which is responsible for the cosmic acceleration today. 



February 1, 2008 7:2 WSPC/INSTRUCTION FILE yykeum 



2 Yong- Yeon Keum 

2. Three possible solutions for Accelerating Universe: 

Recent observations with Supernova la type (SNIa) and CMB radiation have pro- 
vided strong evidence that we Hve now in an accelerating and almost flat universe. In 
general, one believes that the dominance of a dark-energy component with negative 
pressure in the present era is responsible for the universe's accelerated expansion. 
However there are three possible solutions to explain the accelerating universe. The 
Einstein Equation in General Relativity is given by the following form: 

Gf^i, = Rfiiy - -Rgf_i,y = SnGTf^i, + Ag^^, (1) 

Here, G^j/ term contains the information of geometrical structure, the energy- 
momentum tensor keeps the information of matter distributions, and the last 
term is so called the cosmological constant which contain the information of non-zero 
vacuum energy. After solve the Einstein equation, one can drive a simple relation: 

i? 4ttG , „ X A 

— (p + 3,) + -. (2) 

In order to get the accelerating expansion, cither cosmological constant A {lu\ = 
P/p = —1) becomes positive or a new concept of dark-energy with the negative 
pressure (w^ < —1/3) needs to be introduced. Another solution can be given by 
the modification of geometrical structure which can provide a repulsive source of 
gravitational force. In this case, the attractive gravitational force term is dominant 
in early stage of universe, however at later time near the present era, repulsive term 
become important and drives universe to be expanded with an acceleration. Also 
we can consider extra-energy density contributions from bulk space in Brane- World 
scenario models, which can modify the Friedmann equation as oc p + p . In 
summary, we have three different solutions for the accelerating expansion of our 
universe as mentioned in the introduction. Probing for the origin of accelerating 
universe is the most important and challenged problem in high energy physics and 
cosmology now. The detail explanation and many references are in a useful review 
on dark energji^. 

In this paper, we concentrate on the second solution using the quintessence 
field. In present epoch, the potential term becomes important than kinetic term, 
which can easily explain the negative pressure with ~ —1. However there are 
many different versions of quintessence field: K-essencj^, phantorrP^ quintorJl^, 
....etc., and to justify the origin of dark-energy from experimental observations is 
really a difficult job. Present updated value of the equation of states(EoS) are ui = 
— 1.02 ± 0.12 without any supernova 

datjni. 

3. Interacting Dark-Energy with Neutrinos: 

As explained in previous section, it is really difficult to probe the origin of dark- 
energy when the dark-energy doesn't interact with other matters at all. Here we 
investigate the cosmological implication of an idea of the dark-energy interacting 
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with neutrinos 1 ^ I For simplicity, we consider the case that dark-energy and neu- 
trinos are coupled such that the mass of the neutrinos is a function of the scalar field 
which drives th e lat e time accelerated expansion of the universe. In previous works 
by Fardon et alff^l and R. PecceP, kinetic energy term was ignored and potential 
term was treated as a dynamical cosmology constant, which can be applicable for 
the dynamics near present epoch. However the kinetic contributions become im- 
portant to descreibe cosmological perturbations in early stage of universe, which is 
fully considered in our analysis. 



3.1. Cosmological perturbations 

Equations for quintessence scalar field are given by 

^ + 27i0 + a^^^^=O, (3) 



V,s{<P) = V{<i>) + VM , (4) 



Vi(0) - a-^ J ^^q2 + a^ml{cf>)f{q) , (5) 

mv{<P) = iT^iS "p' (as an example), (6) 

where V{4i) is the potential of quintessen ce scal ar field, Vi{(p) is additional potential 
due to the coupling to neutrino particles ^I^^ISl^ and m^((/)) is the mass of neutrino 
coupled to the scalar field. Ti. is d/a, where the dot represents the derivative with 
respect to the conformal time r. 

Energy densities of mass varying neutrino (MVN) and quintessence scalar field 
are described as 



-4 

Pv = a 



3P^ 



















vV+^/o(g) , (7) 

fo{q) , (8) 



(27r)3 



+ V{<j>) , (9) 

(10) 

From equations ([7]) and ([5]), the equation of motion for the background energy 
density of neutrinos is given by 

-f 3n{p, + p.) = ^^r-^'/'(p. - 3P.) . (11) 

0(p 

In our analysis, we are working in the synchronous gauge with line element: 

ds^ = a^ir) [-dr^ + {6ij + hij)dx'dx^] , (12) 

For CMB anisotropics we mainly consider the scalar type perturbations. We intro- 
duce two scalar fields, h{k, r) and 'ri{k, r), in k-space and write the scalar mode of 
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hij as a Fourier integral'- 



hij{x, ''') ^ J d!^f^^ 



(13) 



where k — kk with k^ki — 1. 

The equation of quintessence scalar field is given by 

□0-K=ff(0) =0 . (14) 

Let us write the scalar field as a sum of background value and perturbations around 
it, 4>{x, t) = (('{t) + S(j){x, t). The perturbation equation is then described as 

To describe S (^^^ , we shall write the distribution function of neutrinos with back- 
ground distribution and perturbation around it as 

f{x\T,q,n,)^fo{T,q){l + ^{x\T,q,nj)) . (16) 

After some calculations, we finally obtain the useful equationJI^ 

dVi d In 



d(j) 



■(P.-3P.), (17) 



+^^iSp,-36P^) (18) 
Note that perturbation fluid variables in mass varying neutrinos are given by 



3.2. Boltzmann Equation 

The Boltzmann equation is given in general, 

Dl^dl^dx^dl^dqdl^dn^dl^ fdl 
Dt dr dr dx^ dr dq dr dru \9t 

From the time component of geodesic equation ^1^, 

^P'f = -^^P^P' - ^9°-m,. > (22) 



(21) 



and the relation = a = a "^^/q"^ + a'^mf,, we have 
dq 1,- i 2'rndmdx' 

d7 = -2'^-'^"" «79^^rf^- ^''^ 
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Our analytic formulas in eqs. (|22ll2"3|) are completely different from those of Brook- 
field et alP^j since they have missed the contribution of the varying neutrino mass 
term. In later this term also give an important contribute in the first order perturba- 
tion of the Boltzman equation. The detail calculations will be shown in elsewher j^^l 
The zeroth-order Boltzmann equation is given by 

f=0. (24) 

The Fermi-Dirac distribution 

can be a solution. Here gs is the number of spin degrees of freedom, hp and /cb 
are the Planck and the Boltzmann constants. We assume that MVNs are decoupled 
from the thermal bath when they are extremely relativistic so we can simply replace 
e in the unperturbed Fermi-Dirac distribution by q. Thus we have 

whish can also be a solution of eg. ([24]) . 
The first-order Boltzmann equation is 

Following previous studies, we shall assume that the initial momentum dependence 
is axially symmetric so that ^E" depends on q = only through q and k n. With this 
assumption, we expand the perturbation of distribution function, ^i, in a Legendre 
series, 

*(fc, n, g, r) - ^(-^)^(2^ + l)*,(fc, g, T)P,{k ■ n) . (28) 

Then we obtain the hierarchy for MVN 

<? , T h din fo _ 

= -^fc*i + --— ^ , (29) 
e b omq 

^1 = i^A:(^'o-2*2) + K , (30) 
3 e 

», = _?M2*,-»,)-(-/. + -,)^, (31) 



where 



^ ^ Alk^-^S^^^^^^ . (33) 
3 e q"^ d<f> dlnq 
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Here we used the recursion relation 

{^ + l)p,+l{^i) - {2l+\)^lP,{^x)-m^l{^J) ■ (34) 

We have to solve these equations with a g-grid for every wavenumber k. 
3.3. Quintessence potentials 

To determine the evolution of scalar field which couples to neutrinos, we should 
specify the potential of the scalar field. A variety of quintessence effective potentials 
can be found in the literature. In this paper we examine three type of quintessential 
potentials. First we analyze what is a frequently invoked form for the effective 
potential of the tracker field, i.e., an inverse power law such as originally analyzed 

OA 

by Ratra and Peebles'^, 

V{(t>) = m4+"(/.-" (Model I) , (35) 

where M and a are parameters. 

We will also consider a modified form of V (0) as proposed by Brax and Martin ^31] 
based on the condition that the quintessence fields be part of supergravity models. 
The potential now becomes 

= M'i+"0-"e3'^'/2m^i (Model II) , (36) 

where the exponential correction becomes important near the present time as — *■ 
mpi. The fact that this potential has a minimum for = -^/a/Smpi changes the 
dynamics. It causes the present value of w to e vol ve to a cosmological constant much 
quicker than for the bare power-law potential ^221. In these models the parameter M 
is fixed by the condition that Vi^ « 0.7 at present. 

We will also an alyz e another class of tracking potential, namely, the potential 
of exponential type 

V{<t)) = M'^e"""'' (Model III) , (37) 

This type of potential can lead to accelerating expansion provided that a < ^/2. 
In figure ([1]), we present examples of evolution of energy densities with these three 
types of potentials with vanishing coupling strength to neutrinos. 



3.4. Time evolution of neutrino mass and energy density in scalar 
field 

For an illustration we also plot examples of evolution of energy densities for inter- 
acting case with inverse power law potential (Model I) in Fig. ([2]). In interacting 
dark energy cases, the evolution of the scalar field is determined both by its own po- 
tential and interacting term from neutrinos. When neutrinos are highly relativistic, 
the interaction term can be expressed as 

^^{Py -3P^)~ ;^(an^.')Vi'.„„.ic.. , (38) 



I 



I 
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Fig. 1. Examples of the evolution of energy density in quintessence and the background fields 
as indicated. Model parameters taken to plot this figure are a = 10, 10, 1 for model I, II, III, 
respectively. The other parameters for the dark energy are fixed so that the energy densities in 
three types of dark energy should be the same at present (left-handed side figure). 

Fig. 2. Examples of the evolution of energy density in quintessence and the background fields 
in coupled cases with inverse power law potential (Model I). Model parameters taken to plot this 
figure are « = l,/3 = l,3as indicated. The other parameters for the dark energy are fixed so that 
the energy densities in three types of dark energy should be the same at present(right-hand side 
figure) . 
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Fig. 3. Examples of the time evolution of neutrino mass in power law potential models (Model 
I) with a = I and /3 = (black solid line), /3 = 1 (red dashed line), /3 = 2 (blue dash-dotted line), 
= 3 (dash-dot-dotted line). The larger coupling parameter leads to the larger mass in the early 
universe. 



where Pi^^^^^^i^^^ denotes the energy density of neutrinos with no mass. The term 
roughly scales as cx a^^, and therefore, it dominates deep in the radiation dominated 
era. However, because the motion of the scalar field driven by this interaction term 
is almost suppressed by the friction term, —SHc/). The scalar field satisfies the slow 
roll condition similar to the inflation models, —3Ti.4> ~ ^^^^(Pi^ ~ 3Py). Thus, 
the energy density in scalar field and the mass of neutrinos is frozen there. These 
behaviors are clearly seen in Figs. ^ and ([3]). 
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Fig. 4. (left-hand side figure) The CMB angular power spectra for Model I. The solid line is the 
best fit for the model {{a, (3) = (2.97, 0.170)), the other lines are models with different parameter 
value of a and (3 as indicated. The points are WMAP three year data. 

Fig. 5. (right-hand side figure) The CMB angular power spectra for Model III. The solid line is 
the best fit for the model ((o, /3) = (0.78, 0.28)), the other lines are models with different parameter 
value of a and /3 as indicated. The points are 2dF data. 

3.5. Constrains on the MaVaNu parameters 

As was shown in the previous sections, the couphng between cosmological neutrinos 
and dark energy quintessence could modify the CMB and matter power spectra sig- 
nificantly. It is therefore possible and also important to put constraints on co upling 
parameters from current observations. For this purpose, we use the WMAP3 ^^SEll 
and 2dF-2S.data sets. 



The flux power spectrum of the Lyman-a forest canbe used to measure the 
matter power spectrum at small scales around z < i I34 | 35| has been shown, 
however, that the resultant constraint on neutrino mass can vary significantly from 
X^'Tijy < 0.2eV to 0.4eV depending on the specific Lyman-a analysis used . The 
complication arises because the result suffers from the systematic uncertainty re- 
garding to the model for the intergalactic physical effects, i.e., damping wings, 
ionizing radiation fluctuations, galactic winds, and so on^^. Therefore, we conser- 
vatively omit the Lyman-a forest data from our analysis. 

Because there are many other cosmological parameters than the MaVaNu pa- 
rameters, we follow the Markov Chain Monte Carlo(MCMC) global fit approach 
1221 to explore the likelihood space and marginalize over the nuisance parameters to 
obtain the constraint on parameter (s) we are interested in. Our parameter space 
consists of 

P = {^bh?, VL^h?,H, r, ^, n„ m„ a, (3) , (39) 

where LObh? and Q.ch^ are the baryon and CDM densities in units of critical density, 
H is the hubble parameter, r is the optical depth of Compton scattering to the last 
scattering surface, and are the amplitude and spectral index of primordial 
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Fig. 6. Contours of constant relative probabilities in two dimensional parameter planes for inverse 
power law models. Lines correspond to 68% and 95.4% confidence limits (left-hand side figure). 

Fig. 7. Same as Fig. ||6}, but for exponential type models (right-hand side figure). 

density fluctuations, and {mi,a,(3) are the parameters of MaVaNu defined in sec- 
tions 3.1 and 3.3. We fiave put priors on MaVaNu parameters as a > 0, and /3 > 
for simplicity and saving the computational time. 

Our results are shown in Figs.® - fT]). In these figures we do not observe the 
strong degeneracy between the introduced parameters. This is why one can put 
tight constraints on MaVaNu parameters from observations. For both models we 
consider, larger a leads larger w at present. Therefore large a is not allowed due to 
the same reason that larger w is not allowed from the current observations. 

On the other hand, larger /3 will generally lead larger in the early universe. 
This means that the effect of neutrinos on the density fluctuation of matter becomes 
larger leading to the larger damping of the power at small scales. A complication 
arise because the mass of neutrinos at the transition from the ultra-relativistic 
regime to the non-relativistic one is not a monotonic function of /3 as shown in 
Fig.([3|). Even so, the coupled neutrinos give larger decrement of small scale power, 
and therefore one can limit the coupling parameter from the large scale structure 
data. 

One may wonder why we can get such a tight constraint on /3, because it is 
naively expected that large /3 value should be allowed if il^h^ ^ 0. In fact, a 
goodness of fit is still satisfactory with large /3 value when il^h^ ~ 0, as shown in 
Fig.®. However, the parameters which give us the best goodness of fit does not 
mean the most likely parameters in general. In our parametrization, the accepted 
total volume by MCMC in the parameter space where U^h^ ^ and /3 > 1 was 
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Table 1. Global analysis data within Icr deviation for different types of the quintessence 
potential. 



Quantites 


Model I 


Model II 


Model III 


WMAP-3 data (ACDM) 




2.21 ±0.07 


2.22 ±0.07 


2.21 ±0.07 


2.23 ±0.07 


^CDM h?[W^] 


11.10 ±0.62 


11.10 ±0.65 


11.10 ±0.63 


12.8 ±0.8 


Ho 


65.97 ±3.61 


65.37 ±3.41 


65.61 ±3.26 


72 ±8 




10.87 ±2.58 


10.89 ±2.62 


11.07 ±2.44 




a 


< 2.63 


< 7.78 


< 0.92 







< 0.46 


< 0.47 


< 0.58 




ns 


0.95 ±0.02 


0.95 ±0.02 


0.95 ± 0.02 


0.958 ±0.016 




20.66 ± 1.31 


20.69 ± 1.32 


20.72 ± 1.24 




Oq[102] 


68.54 ±4.81 


67.90 ± 4.47 


68.22 ±4.17 


71.6 ±5.5 


Age/Gyrs 


13.95 ±0.20 


13.97 ±0.19 


13.69 ±0.19 


13.73 ±0.16 




< 0.44 


< 0.48 


< 0.48 


< 1.97(95%C.L.) 


T 


0.08 ± 0.03 


0.08 ± 0.03 


0.09 ± 0.03 


0.089 ± 0.030 







small, meaning that the probability of such a parameter set is low. 

We find no observational signature which favors the coupling between MaVaNu 
and quintessence scalar field, and obtain the upper limit on the coupling parameter 
within 2(7 ranges as 

(3 < 1.11, 1.36, 1.53 , (40) 

and the present mass of neutrinos is also limited to 

^^I'^t^oday < 0.0095, 0.0090, 0.0084 , (41) 

for models I, II and III, respectively. When we apply the relation between the total 
sum of the neutrino masses and their contributions to the energy density of 
the universe: n,yh^ = M,y/ {93.14eV), we obtain the constraint on the total neutrino 
mass: M^, < 0.87eV{95%C.L.) in the neutrino probe dark-energy model. The total 
neutrino mass contributions in the power spectrum is shown in Fig [8l where we 
can see the significant deviation from observation data in the case of large neutrino 
masses. 



In summary, we investigate dynamics of dark energy in mass-varying neutrinos. 
We show and discuss many interesting aspects of the interacting dark-energy with 
neutrinos scenario: (1) To explain the present cosmological observation data, we 
don't need to tune the coupling parameters between neutrinos and quintessence 
field, (2) Even with a inverse power law potential or exponential type potential 
which seem to be ruled out from the observation of to value, we can receive that 
the apparent value of the equation of states can pushed down lesser than -1, (3) As 
a consequence of global fit, the cosmological neutrino mass bound beyond ACDM 
model was first obtained with the value ^ < 0.87 ey(95%CL). 
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10 I ' ' — ^ — ' — -J — ' — ' — ^ ' ' — J 10 I ' ' — ^ — ' — -J — ' — ' — ^ ' ' — J 

1e-05 1e-04 0.001 0.01 0.1 1 1e-05 16-04 0.001 0.01 0.1 1 

k/h [Mpc 'l k/h [Mpc 'l 

Fig. 8. Examples of the total neutrino mass contributions in power spectrum with M„ = 
0.9 ey( left- hand side graph) and with Mi, = 0.3 eV (right-hand side graph). Here the variable 
A is equal to a. 
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